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ABSTRACT 

We analyse a set of new pencil-beam galaxy redshift data in three small regions 
around the South Galactic Pole (SGP) area. We investigate whether we can find any 
evidence of the quasi-periodic peaks discovered by Broadhurst et al. (1990) in the dis- 
tribution of galaxies along the NGP-SGP directions. We use both a power spectrum 
analysis and a cross-correlation with a sliding comb-like window (the comb-template 
technique). Despite the data are less deep (~ 600/i^^Mpc) and certainly not opti- 
mal for such an investigation, there is evidence of the same preferential ~ 130 
Mpc scale in two fields displaced respectively 15° and 30° West of the Broadhurst et 
al. original probe. Taken alone, however, this scale would not be statistically distin- 
guishable from a noise fluctuation. Nevertheless, the statistical significance raises to 
~ 99% when one refers to the conditional probability of finding a peak around the 
same scale measured by Broadhurst et al. . 

Key vifords: Redshift Surveys - Galaxies - Large Scale Structure 
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1 INTRODUCTION 

Deep "pencil-beam" surveys (angular size ~ 1°, depth 
~ 10^ Mpc) (Koo & Kron 1987; Broadhurst et al. 1988) 
have produced important results about the large-scale struc- 
ture of the Universe. They are complementary to the wide- 
angle surveys (100° or larger, but with depth only 100 — 
200 h~^Mpc; e.g. GiovaneUi & Haynes 1991; Strauss & 
Willick 1995; Guzzo 1996). These narrow shots through deep 
space provide a confirmation of strong inhomogeneities in 
the galaxy distribution up to scales of 50 — 100 /i^^Mpc. 

One of the most exciting results obtained from pencil- 
beam surveys has been the detection, along the NGP-SGP 
axis, of a typical scale in the distribution of galaxy con- 
centrations, corresponding to a characteristic separation of 
128 h-^Mpc (Broadhurst et al. 1990, hereafter BEKS). This 
result is consistently confirmed by extensions of the original 
data set (Broadhurst et al. 1995), by independent galaxy 
surveys in nearby areas, as the ESO Slice Project (ESP, 
Vettolani et al. 1995) , and by the large-scale distribution of 
clusters of galaxies (Bahcall 1991; Guzzo et al. 1992; TuUy 



* Based on data collected at the European Soutfiern Observatory, 
La Silla, Cliile. 



et al. 1992). This latter evidence shows how at least the 
first maxima do coincide with large-scale 3D structures, and 
are not artifacts produced by an aliasing effect as claimed 
by Kaiser & Peacock (1990). However, what the periodicity 
precisely means in terms of large-scale structure is still con- 
troversial (e.g. Dekel et al. 1992). A similar effect would be 
expected in only 15% of all directions in phenomenological 
Voronoi models (e.g. van de Weygaert 1992). 

A clearer picture will probably arise only from new 
pencil-beam observations in different directions (see Broad- 
hurst et al. 1995), together with a detailed knowledge of 
the 3D power spectrum of density fluctuations on A > 
100 /i~^Mpc. With this Letter, we try to add a small contri- 
bution in this direction, by analysing some new redshift data 
in the SGP area. All through the paper, we use Ho = 100 h 
km sec~^ Mpc~^ and fio = 1. 



2 DATA GLOBAL PROPERTIES 
2.1 The data 

We use the redshift data discussed in Ettori et al. (1995, Pa- 
per 1 hereafter), where full details about the observing strat- 
egy and the data reduction procedures can be found. These 
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redshifts were collected using the 33-arcmin field ESO fibre 
spectrograph OPTOPUS m throe regions around the SGP, 
named PW, PL and FD, for each of which we summarise 
here the salient aspects. The completeness of each field at 
different magnitude limits is given in Tab. 1. 
Field PW contains 53 redshifts, in the area comprised be- 
tween 22''35'" < a < 22''40'" and -29° < 6 < -27° 45', 



Fig. 1 shows the resulting (j>{dc) for the three regions 
plotted on top of the observed distribution in bins of 30 
Mpc (shaded histogram). In the same figures we show also 
the result of normalizing the data to the total number of 
galaxies expected from the parent photometric catalogue, 
again considering separately the completeness for each bj- 
magnitude bin (open histogram) . 



Field PL contains 415 redshifts, limited to bj < 19.62. 

It extends bectwon 23''39™ < a < 23''52™ and -29° < 
5 < —27° 45', centred on the rich cluster Klemola 44 (Abell 
4038). With a mean velocity < cz >~ 8800 km s"\ the 
cluster overdensity dominates the redshift distribution in 
this beam, contributing ~ 150 galaxies (34% of the total). 
In fact, the main aim of the original survey was to study 
in detail the kinematics and structure of K44 and its sur- 
roundings (Ettori et al. 1996). Despite this bias towards a 
specific redshift peak, the faint magnitude limit of the sur- 
vey provides a large number of foreground and background 
redshifts. To exclude the cluster effect, we also analyse sep- 
arately the subsample of 288 galaxies laying at cz > 11500 
km s~^, that we call PL-K44. 

Field FD is composed by 5 not-contiguous Optopus fields 

(see Paper 1), located between 2''01'" < a < 2^16"" and 
-30° 10' < S < -28° 15', ~ 30° East of the BEKS probe, 
containing 54 galaxies to bj < 19.5. 



2.2 The Selection Function 

The comoving distance for all the galaxies in the sample is 
evaluated as 



dc 



with go = 0.5. The results discussed in the following do not 
change significantly using qo = 0.1. The selection function 
can be then calculated by integrating the Schechter lumi- 
nosity function. 



$(L)d(L) 



(2) 



to the magnitude limit of each field. However, in order to 
take properly into account the magnitude-dependent incom- 
pleteness evident from Tab. 1, we estimate the actual selec- 
tion function separately for four apparent-magnitude bins. 

In other words, wc normalize to a different <j}* within each 
bin, and then calculate the global selection function as 



0(dc) = 



M,+ i(63,do) 



*($*,M)dM. 



(3) 



M,(bj,da) 



For the shape parameters of the luminosity function, 
we adopt the very recent estimate from the ESP survey 
(Zucca et al. 1996), that finds M* = -19.57 and a = -1.22 
. This estimate of the luminosity function is optimal for our 
data, having been performed on same parent photometric 
catalogue (the EDSGC) and to a very similar magnitude 
limit. The selection function also takes into account the K- 
correction, calculated following Shanks et al. (1984), assum- 
ing an Early- type/Late- type ratio of 0.3/0.7. This hypothe- 
sis is necessary, the morphological type being unknown for 
most of our galaxies. 



2.3 Data SpEirseness 

Before starting to look for possible periodic patterns in our 
data, we asked ourselves whether the data were at least 
good enough to answer the zero-order question: "Is the 
data distribution significantly different from a random dis- 
tribution observed with the same selection function?" . One 
way to quantify this difference is to calculate the so-called 
Double- Root Residuals (DRRs). This method compares an 
observed distribution function with a model by estimating 
properly weighed residuals (Gebhardt & Beers 1991). Resid- 
uals are normally calculated giving equal weight to both 
low-populated bins, generally belonging to the tails of the 
distribution, as to central bins containing many more data 
points. Since the residuals depend strictly on ^/N /N fluc- 
tuations, they are not a correct diagnostic of the agreement 
between the data and the model all over the range covered 
by the distribution. The DRR method introduces a square- 
root trasformation, which appropriately weighs the residuals 
of the different bins, showing where there is a significant de- 
viation of the data from the model. If Bi is the data value 
of bin i and Mi is the corresponding model value, the DRRs 
are defined as: 



\qoZ + {qo-\) -l-hV2?o« + l }, (1) I 



DRR^ 
DRR^ 



V2 + ABi - y i AMi if > 1 

1 - VI + 4Mi if = ^ ' 



The local deviations of the DRRs for the given model 
(assumed to describe sufficiently well the global trend of 
the data), give directly a confidence level in terms of stan- 
dard deviations. To reject the null hypothesis at the 95% 
confidence level, for example, we will require bins with 
|DRR| > 2. 

Fig. 2 shows the DRR's for our data, when compared to 
the selection function. For PW and FD the null hypothesis 
can be rejected at the 95% confidence level only for a few 
bins. This already tells us not to expect to obtain much 
significant statistical information from these two data sets. 
On the contrary, PL shows rather significant fluctuations, 
with peaks in corrispondence to K44 and at 170 and 570 
ir^ Mpc and troughs around 250, 350 and 450 ir^ Mpc. 



3 PERIODICITY ANALYSIS 

Having established that the data show at least some devi- 
ations from an uniform distribution, we address now the 
question of the presence within the three regions of any 
significant periodic pattern similar to that originally found 
by BEKS. To this end, we apply two statistics, the clas- 
sical Power Spectrum analysis, and the so-called Comb- 
Template Test. 
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3.1 The Power Spectrum analysis 



3.2 The Comb-Template Test 



Considering a set of N data with comoving distance dj , 
define the power spectrum as the quantity 



Pk = 



(5) 



This can bo seen as the squared amplitude of a discrete 
Fourier transform of Soirac functions, each located at 
position dj corresponding to the galaxy positions. Under the 
null hypothesis of a uniform galaxy distribution on scales > 
30 h~ Mpc (Szalay et al. 1991), the probability distribution 
of Pk is equal to a x^(2)/2 distribution, i.e. an exponential 
(since the real and imaginary parts of the Fourier coefficients 
are Gaussian distributed and statistically independent; cf. 
Burbidge & O'Dell 1972; see Amendola 1994 for the non- 
Gaussian case). The expectation value of this distribution 
is unity. To explore the significance level of any power Pk 
present in the spectrum, let us consider the new variable 
Zk = Pk/ < Pk >, where < > is the mean power in 
the spectrum. The observed distribution of Zk will have a 
mean which is 1 by defirntion. We can then evaluate the 
significance of any observed power by comparing it to the 
null hypothesis, still described by an exponential. 

It can be shown (cf. Burbidge & O'Dell 1972, Szalay 
et al. 1991) that for K independent test wavenumbers the 
significance level of any peak in the power spectrum is 
given by 

p {zk > zl) = 1 - [1 - exp(-4)] . (6) 

To establish the significance of a peak, we need therefore to 
estimate the two parameters < Pk > and K for each field. 
< Pfc > is the mean power estimated on the K wavenumbers; 
K can be calculated as 



K ■ 



Ak 



(7) 



where k,nax and kmin correspond respectively to the inverse 
of the minimum and maximum scale that we chose to anal- 
yse (see below). Ak is the spacing in wavenumber beetwen 
two contiguous fci's, i.e. the smallest independent frequency 
and/or the largest comoving range, that for our case cor- 
responds to Ak = J- — — = jj^. For our data, we 

limit the range of investigation to scales between 30 and 
400 Mpc, i.e. within a range where we can neglect the 
effects of small-scale clustering ( "red noise" ) and of the un- 
certainty on the selection function at large distances. These 
yield k^ax = 0.0333 /iMpc"^ and fc^™ = 0.0025 /iMpc"^ 

To produce a smoother picture of the power spectrum, 
we oversampled the number of wavenumbers, as shown in 
fig. 3, i.e. calculating the power at "dependent" wavermrri- 
bers separated by less than Afe — 1/700. It should be clear 
that this oversampling does not add any contribution to the 
significance level of an investigated power. 

The final results are shown in fig. 3 and listed in Tab. 2. 
None of the peaks observed in the different fields provides 
a statistically significant evidence for a preferred clustering 
scale. The probability that the highest peaJi is simply a ran- 
dom fluctuation of the noise is everywhere larger than 20%. 
However, it is worth noticing how for PL-K44 and PW the 
value of the most probable period is interestingly close to 
the BEKS value of 128 /i"^ Mpc. 



The Comb- Template Test (CT hereafter), has been used by 
Duari et al. (1992) to look for periodicity in the distribu- 
tion of QSOs. Its main advantage is to be able to identify, 
in a simple way, the starting point of a periodicity (i.e. its 
phase), providing also a relative probability value for the 
more significant period. 

Assuming that the objects are concentrated within 
peaks of width w, spaced with a period D, we slide a comb- 
like template with periodic "teeth" across the galaxy redshift 
distribution and evaluate the function 



c{di, da) 



jl ff 
\ o1 



di - {do + D nint {^^^)}\ < f 
otherwise 



(8) 



where do is the starting point and nint{x) is the nearest 
integer to x. 

For a discrete distribution of doirac functions, the test 
consists in evaluating the correlation function 



(9) 



For an uniform galaxy distribution, it can be shown that 
t{do) has the a mean value fj, = Nw/D and standard devi- 
ation at = \/ n{i — w/D). The highest value of t{do), tmax, 
corresponds to the most probable starting point do. In ad- 
dition, given a test period the null hypothesis of a uniform 
distribution can be rejected at the n-a confidence level if 
tmax is larger than n by the same amount. Fig. 4 shows the 
behaviour of t{do) for our 4 samples, when the D parame- 
ter is fixed to the value obtained from the power spectrum 
analysis and the width w is set to 15 Mpc. We did not find 
any significant dependence of the results on variations of 
this latter parameter. Note that the do is measured from 
the lower limit of the data set in comoving distances, which 
for the PL-K44 sample corresponds to the redshift cut at 
cz = 11500 km s-^ 

A second interesting step is to compare the values of 
tmax obtained for different periods, to find which period cor- 
responds to the "least random" distribution. The results of 
this test, performed on scale between 50 and 200 Mpc, are 
shown in fig. 5. This analysis is conceptually different from 
the power spectrum, but, as can be seen in the figure, it 
provides a substantial confirmation of the results of the har- 
monic analysis. Only in one case (PL-K44), there is some 
evidence for a different specific scale at 73 ~ 70 Mpc, which, 
interestingly, has a similar counterpart also in field FD. We 
summarise the results of the CT test in tab. 3. 



4 DISCUSSION AND CONCLUSIONS 

Let us discuss first the case of the field PL, where the situa- 
tion is complicated by the dominating overdensity of the K44 
cluster of galaxies. In fact, the power spectrum (fig. 3, upper 
panels) shows a rmmber of peaks, with no particular pref- 
erence for the BEKS scale. The main "period" which arises 
from both the power spectrum and the CT test (fig. 5), 
clearly corresponds to the separation between the cluster 
and the second redshift peak, ~ 98 /i~^Mpc. The critical 
effect of the cluster on the harmonic analysis is also evi- 
denced by the value of the starting point provided by the 
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CT test - i.e. the phase of the hypothetical "periodicity" 
- which corresponds exactly to the position of the cluster 
(85 — 86 h~^Mpc). It is interesting to note, however, that 
when the cluster is removed from the sample, a peak arises 
in Pk right at the BEKS scale, although it is just slightly 
more than a ~1-(t fluctuation and therefore not statistically 
significant. 

The field PW is the only one with a dominant peak 
in the power spectrum, corresponding to a period ~ 
120 h~^Mpc. Also in this case, however, its significance 
is low, due to the sparseness of the data evidenced by the 
DRRs: the null hypothesis can be rejected only at the 77% 
confidence level. The CT test shows the same tendency for 
a preferential scale around a similar value. 

In the field FD there is some consistency for a preferred 
scale at ~ 70 h~^Mpc from both the Pk and CT analyses. 
However, the significance of this peak is even lower than for 
the other fields (with a 54% confidence level). We mention 
it for two reasons: first, because the CT test shows a similar 
scale also for the field PL-K44; second, because a similar 
scale was detected by Mo et al. (1992) in their analysis of a 
number of different samples, and is close to a kind of second 
harmonic of the BEKS scale. 

So far we have considered the simple absolute signifi- 
cance of the potential periods evidenced by our analysis. In 
other words, we have asked ourselves whether in our three 
sets of data there is any significant evidence for a regularity 
with a preferred scale in the galaxy one-dimensional dis- 
tribution. The answer to this question is clearly negative: 
we have shown how our set of data is essentially not good 
enough on its own for this purpose. 

However, there is a more specific question that we can 
ask, i.e. given the three independent data sets, is there any 
evidence for a peak in the power spectrum around the same 
value found by BEKS, and how significant is it? The answer 
to this second question is more interesting, as we have seen 
that in two directions, PL and PW, there is indeed a peak 
around that same scale. The significance level of this event 
will be now given by the conditional probability of the two 
events of (1) having a peak of that height anywhere in the 
explored range, and (2) finding it around 130 /i~^Mpc. In 
practice, this results is the product of the absolute proba- 
bility estimated in table 2, with the inverse of the number 
of frequencies explored, which is 22 in our case. This gives 
a probabUity of 1.4% for PL and 1% for PW (almost 99% 
confidence level) that a scale of 120-130 h~^Mpc is shown 
in a sample of galaxies randomly drawn from a population 
of uniform distribution. 

It seems difficult to ascribe the systematic detection of 
this same scale from many different samples, including the 
rather sparse one examined here, to chance coincidences. 
While we were writing this paper, evidence for excess power 
on 100 h^^Mpc scale was also found in a 2D power spec- 
trum analysis of the Las Campanas Redshift Survey (Landy 
et al. 1996), further supporting the idea that wavenumbers 
in the ~ 100 — 200 /i~^Mpc range have a particular meaning 
in the structuring of our Universe. One natural possibility 
is simply that both ID (skewers) and 2D (slices) projec- 
tions amplify and enhance a true maximum in the three- 
dimensional power spectrum. This maximum is expected in 
all variants of CDM models, as the signature of the horizon 
size at the epoch of matter-radiation equivalence. Observa- 



tionally, the spectrum has to turnover around these scales, 
to be able to connect the very large-scale estimates from mi- 
crowave anisotropy measurements (e.g. Hu et al. 1996), to 
the ~ 100 /!,~^Mpc-scale direct measurements from galaxy 
redshift surveys (e.g. Lin et al. 1996). To produce the ob- 
served quasi-regular distribution of objects, the spectrum 
turnover has to be rather sharp (Frisch et al. 1995). Intu- 
itively, the sharper the 3D peak is, the better a well defined 
preferential cell size is defined. The next generation of red- 
shift surveys of galaxies and clusters (see e.g. Guzzo 1996 for 
a review), will provide a first possibility to study both ID 
skewers and 3D power spectra from the same survey volume, 
thus allowing to finally clarify the issues discussed here. 
Acknowledgments. We thank G. Chincarini for his com- 
ments suggestions during the early phases of this work. SE 
thanks A. Gaspani for useful discussions on power spectra. 
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Figure Captions 

Figure 1. a) Redshift distribution in 30 Mpc bins 

in tlic PL field. Tiic siiadcd liistogram is the observed dis- 
tribution, while the open histogram shows the result of nor- 
malizing each apparent magnitude bin for the completeness, 
as shown in Table 1 and discussed in the text. The solid 
curve is the selection function calculated using eq. 3. Verti- 
cal dashed lines indicate the best periodicity as obtained in 
chapter 3. b) The same, for field PW. c) The same, for field 
FD. 

Figure 2. The results of the DRRs analysis, showing 
the significance of the fluctuations in the data with respect 
to the null hypothesis of a uniform distribution with the 
same selection function. The dashed lines at ± 2 indicate 

the 95% significance level. 

Figure 3. Power Spectra for the four data sets anal- 
ysed. The vertical dashed lines correspond to the 128 
Mpc BEKS periodicity (fc ~ 0.0078 h Mpc'^). 

Figure 4. The distribution of the phases do from the 
Comb-Template Test. The period adopted (as obtained 
from the power spectrum), is reported for each data set. 

Figure 5. The distribution of the test periods from the 
Comb-Template Test, n gives the number of a deviations 
from a uniform distribution, for any given D. The arrows 
mark the highest peak, corresponding to the best period. 
The BEKS periodicity is given by the dashed lines. 
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Table 1: Completeness in the three fields at different magnitude Umits. Ustars 
is the nmiiber of objects found to be misclassified stars in the EDSGC and 
therefore excluded a posteriori from the sample (see Paper 1). is the actual 
number of galaxy redshifts available. Note that the completeness is estimated 
in a conservative manner, i.e. assuming that the remaining unobserved objects 
in the target catalogue are all galaxies. If one applied the stellar contamination 
percentages to the entire parent catalogues, the observed completenesses would 
improve of a few percents, essentially in the fainter magnitude bins. 
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Table 2: Results of the power spectrum analysis. Tpeak is the period corre- 
sponding to the highest peak detected. Ppeak is the height of the peak. < Pk > 
is the mean noise level, and p is the probability that the peak is a random noise 
fluctuation. 
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3.87 


71 


4.00 



Table 3: Results of the Comb-Template Test, do is the most probable value 
for the phase (i.e. the starting point), of a periodic pattern having the period 
as yielded by the power spectrum. Ua- is the confidence level (in number of 
(t's), with which the null hypothesis of a uniform distribution can be rejected. 
Dmax is the most probable period obtained by the test when do is left as a free 
parameter. Ua-^max is its confidence level, again in terms of standard deviations, 
as defined in the text. 




20 40 60 80 100 120 20 40 60 80 100 120 140 




20 40 60 20 40 60 80 100 120 



(Mpc) 



